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Recent high-sensitivity observation of the nearby radio galaxy M87 have provided im- 
i ^ | portant insights into the central engine that drives the large-scale outflows seen in radio, 

optical and X-rays. This review summarizes the observational status achieved in the 
high energy (HE;< 100 GeV) and very high energy (VHE; >100 GeV) gamma-ray do- 
mains, and discusses the theoretical progress in understanding the physical origin of this 
emission and its relation to the activity of the central black hole. 
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£\J 1. Introduction 

The advanced capabilities of current gamma-ray instruments offer a unique tool to 
, ^ elucidate the nature of the central engine in extragalactic objects. At the time of 

writing, about 1000 active galaxies have been detected above 100 MeV by Fermi- 
LAT. 2 Ground-based gamma-ray instruments on the other hand have established 
TeV emission and spectra for about 50 active galaxiesj^] The majority of these are 
blazars, only a minor fraction are radio galaxies. The giant elliptical Virgo-cluster 
galaxy M87 (NGC 4486, 3C274) is particularly unique among those sources, being 
the nearest AGN in the northern sky, one of the first extragalactic objects dis- 
covered (in 1918) to have an optical jelP, the first extragalactic X-ray source to 
be identified^, and the first non-blazar active galaxy detected at VHE energies.^ 
Because of its proximity, many astrophysical phenomena can be studied in excep- 
tional detail, which has led M87 to become a prominent benchmark for theoretical 

a For an up-to-date list see the online TeV Source Catalog (TeVCat) at |http:/ /tevcat.uchicago.edu/| 
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research. In particular, since the jet in M87 is substantially inclined, it is fairly easy 
to study its core - the vicinity of the black hole. 

2. General properties of M87 

(1) Located at a distance of d ~ 16.7 Mpc (redshift of z = 0.004; 1 mas=0.081 
pc)J33 M87 is commonly counted among the FR I clasiP of (weak power) radio 
galaxies. It is considered to host one of the most massive black holes (BH) in the 
Universe with mass estimates based on gas or stellar dynamics in the range of 
Mbh ^ (2-6) x 10 9 M ) 9 l lo l n l(see also Ref.DJfor a fundamental plane inference) P 
There is evidence in HST data for a possible offset between the galaxy (bulge) photo- 
center and the AGN that might be related to a gravitational wave recoil resulting 
from the coalescence of a supermassive BH binaryP 

(2) M87 is well-known for a bright one-sided, quasi-straight jet (of projected length 
~ 2 kpc) that extends from the center to the outer parts of the source and can be 
clearly seen in radio, optical and X-rays, cf. Fig. ll 15 l 16 l On very large scales, there 
is evidence (VLA, 90 cm) for two directed, coherent outflows powering the giant 
radio halo extending about 40 kpc from the nucleusP 

(3) High-resolution VLBA images suggest that the innermost part of the jet seen at 
radio frequencies is formed very close to the black hole, with an apparent (initial) 
'full opening angle' of approximately 9j <a — 60°. Significant collimation of the jet 
already occurs within a few tens of Schwarzschild radii (r s ) from the black hole, see 
e.g. Fig. 2, but seems not complete until about 1000r s ll- 8 l 19 l 20 l 21 l An apparent full 
opening angle of 60° (i.e., intrinsic full opening angle 0j ~ 9j >a sin 9 ~ 30° for a jet 
inclination 9 = 20°, cf. Ref. [2"2"|) up to projected distances of about 100r s J22l suggests 
that the (radio) jet already has a comparable width on those scales. 

(4) Despite its huge black hole mass, M87 is not a powerful source of radiation. 
Common estimates for its total nuclear (disk and jet) bolometric luminosity do not 
exceed Lboi — 10 42 erg/s by muchU 7 -! 24 ! 26 ! When this is compared with the canon- 
ical Eddington luminosity value ^Edd, M87 becomes highly under-luminous by a 
factor l e = Lboi/^Edd < 3 x 10~ 6 . This has led to the proposal that M87 is a 
prototype galaxy, where accretion occurs in a non-standard, advective-dominated 
(ADAF) mode characterized by an intrinsically low radiative efficiency! 27 ! 28 ! 29 ! In 
such a case, higher accretion rates are permissible, so that for M87 accretion rates 
m up to some ~ 10~ 4 mEdd (in terms of the Eddington rate) become possible^ 
values that are indeed compatible with the (Bondi) upper limit derived from Chan- 
dra observations.^ Assuming equipartition with the accreted matter, characteristic 
ADAF magnetic fields in the vicinity (~ r s ) of the black hole (and thus likely the 
maximum magnetic field feeding the jet) would then be of order B ~ 200 GP^ 

(5) M87 has been a prime target for MHD jet simulation studies (see e.g., Refs. I3"2T 
[37]). In principle, the maximum possible kinetic power Lj carried by a jet is expected 
to be comparable to the source's accretion power ~ tiic 2 (e.g., Ref. 38). For M87 
this translates into an upper limit Lj ~ 3 x 10 44 (m/10~ 4 mEdd) erg/s, which seems 
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consistent with other estimates for the time-averaged bulk kinetic power as inferred 
from its large-scale radio and X-ray morphology, ranging from Lj ~ 3 x 10 up 
to a few xlO 44 erg/sllll 25 ! 39 ! 40 ! Obviously, if accretion would occur in the standard 
(cooling-dominated) mode, current accretion rates would be too small to account 
for the mean kinetic jet power. 

(6) Optical (HST) observations of M87 reveal apparent superluminal motion of jet 
components about 0.5 kpc away from the central black hole, suggesting bulk flow 
Lorentz factors Tj ~ 6 and a jet orientation within 9 ~ 19° to the line of sightP-D 
Radio (20cm VLBA) observations also indicate apparent superluminal speeds up to 
0s ~ 4.3 ± 0.7 of features moving through knot HST-1 (located 0.85 arcsec ~ 60 
pc in projection away)P^l Slower speeds between (0.6 — 0.8)c within 80 pc of the 
nucleus have been reported based on HST dataP^ On the other hand on scales 
of one parsec (~ 1.6 pc) no significant relativistic motion is evident from radio 
(2cm) observations)-^ and on sub-parsec scale only mildly relativistic (radio jet) 
bulk speeds (/? ~ 0.7) and somewhat larger jet inclination angles (yet with velocity 
measurements based on one pair of observations only) have been inferred based on 
43 GHz radio VLBI observations.^ A likely range for the jet inclination thus seems 
to be 9 ~ (15 — 25)°P^1 This would be largely consistent with M87 being a "non- 
blazar" source (or "misaligned BL Lacertae object"^) in the sense that most of its 
jet emission (at least at radio frequencies) is characterized by rather weak Doppler 
factors D = l/[Tj(l — f3cos9)) ~ 3 and moderate Doppler beaming onlyEMl 
Quite surprisingly, given this seemingly unfavorable condition, M87 has turned out 
to be a highly variable TeV gamma-ray emitter, extending the class of VHE emit- 
ting AGN from the classical blazars to radio galaxies. Core features at gamma-ray 
energies are summarized in the following. 

3. Observational findings towards high energies 
3.1. Very high energy (VHE) results 

(1) For the first time, a significant (> 4<r) excess of photons from M87 above energies 
of ~ 700 GeV has been measured by the HEGRA stereoscopic systemP Subsequent 
observations by the current generation of VHE instruments have confirmed gamma- 
ray emission at TeV energies and further demonstrated the exceptional behavior of 
M«7 l 47 ' 48 ' 49 ' 

(2) H.E.S.S. has detected fast (day-scale) TeV flux variability during a high VHE 
source state in March/ April 2005 (average photon flux a factor ~ 5 higher than in 
2004), and demonstrated that the VHE flux of M87 (between 2003-2006) also varies 
on long (year- type) time scalesP^The measured VHE spectrum above 0.35 TeV was 
found to be consistent with a (single) hard power-law (photon index T = 2.22±0.15, 
compared to 2.62 ±0.35 for the 2004 low state), and to extend beyond 10 TeV. Sim- 
ilar to other wavelengths, the TeV output was relatively moderate, with an average 
(for the years 2003 to 2006) isotropic TeV luminosity LtoV — 3 x 10 40 erg/s (the 
luminosity being roughly three times larger for the 2005 high state). 
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(3) VERITAS observations during a rather low state in spring 2007 confirmed the 
hard VHE spectrum and provided further evidence for year-scale variability,^ while 
MAGIC observations^ during another VHE high state in February 2008 again re- 
vealed substantial day-scale variability (cf. also the 4d-VHE flare observed by VER- 
ITAS in 02/2008, Ref.EU). Interestingly, the MAGIC findings seemed to indicate a 
different variability behavior; above 350 GeV considerable flux variability on time 
scales as short as ~ 1 day has been found, while at lower energies (150 — 350 GeV) 
no significant variability was detected.^ This could perhaps point to the appear- 
ance of a new (additional) emission component at the highest energies. 

(4) Day-scale VHE activity with a well-defined rise/fall structure of the flux was 
also observed during another VHE high state (exceeding 10% of the Crab above 
350 GeV) in April 2010. 1 52 VHE flux (exponential) rise and decay times of - 3 d 
and ~ 1 d, respectively, have been reported by VERITAS. See also Fig. [3] 

(5) There are indications for spectral variability in this and the previous flares, 
with photon indices (assuming a single power law) in the low (quiescent) and high 
states closer to T ~ 2.6 and T ~ 2.2, respectively! 47 ! 49 ! 52 ! 

(6) Both the hard VHE spectrum extending beyond 10 TeV, and the observed rapid 
VHE variability down to (flux doubling) time scales of Ai Q b s — 1 day are remarkable 
features for a source whose overall jet is commonly believed to be substantially in- 
clined to the line of sight. The reported VHE flux variations are the fastest observed 
in any waveband from M87 so far, and already imply a size of the 7-ray emitting 
region (R ~ DAi b s c), that is of the same order of magnitude as the Schwarzschild 
radius of the black hole in M87. Given current evidence, this may possibly point to 
the black hole vicinity as the most likely origin of the VHE radiation. 

(7) The restricted angular resolution of current Cherenkov telescopes (~ 0.1°, cor- 
responding to 30 kpc at the distance of M87) however, does not allow to resolve the 
jet in M87, so that the true location of the VHE region in charge usually remains 
ambiguous when based on the VHE results alone. On the other hand, because of 
weak Doppler boosting, there is not much uncertainty concerning the upper limit 
on the size of the VHE source. In addition, for the 2008 (February) VHE high state 
a combined radio (VLBA, 43 GHz, high-resolution 0.21 x 0.34 mas) and VHE mon- 
itoring effort found that the VHE gamma-ray outburst was followed by a strong 
increase of the radio flux very close to the core (within ~ 200r s along the jet for 
a viewing angle 15°; note that there is reasonably experimental evidence that the 
radio core is located close to the black hole, Ref. [5T]). As the flare progressed, the 
region around the core brightened up substantially (associated average /3 8 = 1.1), 
eventually reaching radio flux levels beyond any seen before, suggesting that a new 
component was ejected from the core. If the increase of the radio flux is caused by 
a decrease in synchrotron self-absorption, then efficient particle acceleration must 
take place very close to the black holeP^Note that a similar radio/TeV connection 
has not been found for the 2010 VHE flare. 
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3.2. High energy (HE) results 

At high energies (> 100 MeV), Fermi-LAT reported the detection of gamma-rays 
from M87 up to 30 GeV based on 10 months of data taken between 08/2008- 
05/2009P Given its angular resolution (0.8° at 1 GeV) the HE location is difficult 
to pin down. No evidence for significant HE flux variations over the period of ob- 
servations has been found, the reported HE light curve (lOd-bins) being consistent 
with no variability. It should be noted, however, that given the low HE flux of M87 
and the limited sensitivity of Fermi-LAT, a potential variability on smaller (day- 
scale) timescales usually cannot be probed. 

The corresponding isotropic luminosity in the HE band is L{> 100 MeV) ~ 5 x 10 41 
erg/s. Interestingly, the HE spectrum essentially appears consistent with a power- 
law whose photon index (r = 2.26 ± 0.13) is comparable to the VHE one in the 
high states. However, given the observed HE flux, a simple extrapolation of the 
HE power-law to the VHE regime is not sufficient to account for the flux measured 
during the TeV high states (which, however, only last for a couple of days and seem 
to occur with duty cycles < 10%). Similar results have also been obtained more 
recently based on an analysis of a larger (two year, 08/2008-08/2010) data setP 
Comparing the HE fluxes up to 05/2009 with those measured afterwards, no evi- 
dence for variability in the fluxes above 100 MeV has been found, but there seems 
to be some weak indications for a change in the fluxes above 1 GeV. 1 

3.3. X-ray results 

(1) M87 has been a target for several X-ray telescopes, including the Einstein 
ObservatoryJUEH ROSAT^O an d XMM-NewtorP 1 . However, although the spatial 
resolutions of these instruments were sufficient to detect the jet, observations with 
the Chandra satellite were needed to resolve the jet into more than one knot, and 
to allow a better localization of the X-ray emission. These observations were started 
in July 2000)151^1 with a detailed M87 monitoring program implemented in 2002.^ 
Given Chandra's subarcsec resolution (~ 0.8 arcsec), X-ray emission from the core 
(nucleus) and the innermost knot HST-1 can be distinguished. 

(2) The coincidence of VHE 7-ray flaring activity in 2005 with a giant X-ray flare 
from HST-1 (with isotropic X-ray luminosities up to a several 10 41 erg/s) initially 
prompted speculations that the observed VHE emission may have originated in 
HST-1 as wellPThe fact, that the Chandra X-ray flux from HST-1 did not change 
much during subsequent VHE flaring episodes (lasting from a few days up to a weak 
or two) has been often taken to disfavor HST-1 as a site of the rapid TeV flaring 
activityl 1 l 48 l 71 llt should be noted, however, that from a physical point of view, even 
a possible anti-correlation is, in itself (!), not so much telling, as in a simple inverse 
Compton model a variation of the magnetic field could lead to anti-correlation. 

(3) While Chandra observations during VHE flares have not yet been achieved, the 
nucleus usually shows comparably high X-ray fluxes in Chandra pointing close to 
the VHE flaring states 
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(4) Observed characteristic X-ray variability (flux doubling) time scales for HST-1 
are of the order t v ~ 0.14 yr (~ 50 d)) 61 l 68 l while the core (nucleus) seems to display 
faster variability (t v . c ~ 20 d); there is some evidence that the shortest nuclear (dou- 
bling) time scale that one can measure with Chandra may even be as small as a few 
dayslil^syHI jf t ru6j this would most likely argue against an ADAF (bremsstrahlung) 
origin of the Chandra X-ray variations.^ 

4. Theoretical interpretations 

A variety of models have been proposed in the literature to account for the HE and 
VHE 7-ray emission from M87, with production sites ranging from large (hundred 
pc and beyond; e.g., knot HST-1 or knot A) via the innermost (sub-pc) jet scales 
down to the immediate environment (magnetosphere) of the supermassive black 
hole. Some representative candidate sites and mechanisms will be discussed in the 
following. 

4.1. HST-1 

The production of TeV photons by inverse Compton (IC) scattering requires TeV 
electrons. As these electrons also emit optical to X-ray synchrotron emission, the 
bright optical/X-ray features in M87 are putative sites for the origin of the observed 
TeV emission. A priori, this includes the inner jet (core), the HST-1 complex (at a 
de-projected distance z > 120 pc) and knot A (at a distance of about 1 kpc). How- 
ever, knot A, with a typical resolved X-ray size ~ 80 pcpl is commonly considered 
to be disfavored (based on light travel arguments^! as a significant contributor to 
the observed variable VHE emission on daily and yearly time scales. HST-1 (typical 
radio size ~ 0.2 pc)J^ on the other hand, has been thought to remain a promis- 
ing site where IC up-scattering of ambient starlight photons by shock-accelerated 
electrons could lead to variable TeV emission! 42 ! 65 ! This complex has shown inter- 
esting activities during the past years, including phases of strong optical and X-ray 
brightening and the appearance of superluminal sub-structures H3HEH 
The fact that there seems to be no longterm correlation between the Chandra HST- 
1 X-ray and the TeV flux does not necessarily (as argued above) provide cogent 
evidence against such a scenario. Instead, the observed rapid TeV variability on 
time scale of ~ 1 day, implying a source size of R ~ DcAt - 8 x 10 15 (Ai/l d)(D/3) 
cm, appears much more challenging to account for. In the case of a conical jet 
structure, for example, an unlikely small opening angle 9j = R/z ~ 2 x 10 _5 rad 
would be required.^ Perhaps more likely is a scenario where HST-1 is the site of 
a hydrodynamic re-collimation shock (e.g., Ref. H2J but see also Ref. IMjP^! where 
convergence occurs and part of the jet kinetic energy is converted into radiative 
energy (e.g., X-rays and possibly 7- rays). In order to account for the observed fast 
(days-scale) TeV variability, strong collimation would be needed. Recent calcula- 
tions suggest, however, that due to its low radiative efficiency, only modest colli- 
mation is achievable. While sufficient to account for the X-ray variability observed 
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from HST-1 (t v ~ 0.14 yr), it would not be sufficient to account for day-scale TeV 
variability.^ Moreover, the radio to X-ray spectrum of HST-1 (in 2005) can be suc- 
cessfully modeled as a broad-band synchrotron spectrum with a cooling break at 
Vb ~ 10 15 Hzl 68 l 103 l Plausible value for the magnetic field strength at the position of 
HST-1 are of the order 5 x 10~ 3 G and below, suggesting a characteristic synchrotron 
cooling time scale for electrons up-scattering starlight photons to the TeV regime of 
t' s ~ (10 6 /7')(5 x 10~ 3 /B') 2 yr. This value is in fact close to the observed variability 
time scale in the UVP^It seems hard to imagine that any supposed short-term TeV 
variability of HST-1 will not be smeared out on comparable time scales. Based on 
current evidence, it thus seems that (at least) for the episodes of rapid TeV flaring 
activity an origin in the sub-parsec scale jet and below is preferred. It could well be, 
however, that HST-1, the kpc-scale jetP^ and/or perhaps even knot contribute 
to the overall, quiescent VHE emission in M87. 

4.2. Leptonic (inner) jet models 

In this approach, the (non-simultaneous) SED representation is often interpreted 
as showing a break between 10 12 — 10 13 Hz, and a possible synchrotron cut-off 
around or below 10 18 Hz. If this is indeed true, however, then homogeneous (one- 
zone) synchrotron self-Compton models are (under the proviso of moderate Doppler 
boosting) at best only marginally able to match the quiescent TeV flux in M87 and 
generally unable to account for its flaring statelMMH] 

A number of different alternatives have thus been considered that could allow accom- 
modating the TeV results. These include: Dropping the single zone approach by as- 
suming (1) a longitudinal jet stratification (as in the decelerating jet flow/upstream 
Compton [UC] scattering scenario, Refs. [72J ITS]) . (2) a radial jet stratification (as 
in the spine-shear approach )P2 or (3) multiple emission zones with changing line-of- 
sight (LOS) directions (differential Doppler boosting) | 80 l 83 l 

In models of the type (1) pi the TeV emission turns out to be dominated by inverse 
Compton emission from the fastest (r.,- up to a few tens) part of the jet on tar- 
get (synchrotron) photons mainly produced in the downstream (further away and 
slower, Tj ~ a few) jet parts. While such models bear the potential to resolve some 
BL Lac/FR I unification problems, they are, when applied to M87, not without 
their challenges: The resultant VHE spectra tend to be too soft; there is as yet no 
direct observational evidence for a decelerating flow, in fact, the jet in M87 rather 
seems to be accelerating (see 2. (6) above); moreover, the small Doppler factor (e.g., 
D = 1.5 for Tj — 10 and 9 = 20°) associated with the fastest part of the jet and the 
anticipated scale of the jet base usually makes it quite challenging to accommodate 
the observed short time (day-scale) TeV variability. 

Models of the type (2), on the other hand, rely on a possible radiative interplay 
between a fast inner core (spine) and a slower outer (shear) layer. The underlying 
assumption of a velocity gradient across the (inner) jet is indeed compatible with 
radio VLBA imaged^ showing a limb-brightened structure (cf. also Refs. 1351 [75]) 
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and seems supported by relativistic jet simulations!^ In principle, such a setup 
is known to be favorable to further shear- type particle acceleration.^! Yet, even 
when this is neglected, each component should see the (external) radiation of the 
other boosted, thereby enhancing the inverse-Compton emission of both the spine 
and the layer. Given the supposed inclination of the jet in M87, its infrared-optical 
spectrum is expected to be dominated by synchrotron emission from the weakly- 
beamed (!) spine and its TeV part by (external) inverse Compton emission from the 
stronger-beamed layer .ESI The scenario requires the zones to be (quasi) co-spatial 
for the radiative interplay to work efficiently. This, however, could also lead to some 
problems concerning the observed TeV emission in that a simultaneous fit to the 
broadband spectrum seems to imply a rather dense (weakly de-beamed) infrared- 
optical radiation field resulting in considerable 77-absorption and a steep slope at 
TeV energies. Existing model calculations in fact appear unable, for example, to 
reproduce the hard (2005) TeV spectrum seen by H.E.S.S., and it seems not yet 
obvious how this might best be recovered. 

In ("mini-blob") models of the type (3), small sub-regions (e.g., plasma blobs) are 
taken to follow different trajectories within the jet formation zone (~ 50r s ), thereby 
allowing for situations where the velocity vector of a region becomes close to the 
line of sight, so that for this region strong Doppler boosting effects become possi- 
ble (despite the overall misalignment). Such a scenario could in principle account 
for VHE high states that are accompanied by a dramatic change in X-ray flux (as 
related to a SSC-emitting blob moving along the LOS). Again, the success is not 
without its challenge. Modeling of the VHE high state in 2005,^ for example, in- 
dicates a significant deviation from equipartition, with the pressure in relativistic 
electrons exceeding the one in magnetic field by a factor 3> 10 3 . Moreover, required 
magnetic field strengths in the emitting region (B ~ 10 mG) would be a factor 
~ 100 lower than the value anticipated for the local (jet) field based on magnetic 
flux conservation. Perhaps, some of the restrictions could be relaxed if these blobs 
are also considered to carry a significant azimuthal velocity component as expected 
for a rotating jet flow near its baseJ^il^HI or if the blobs (also called "mini-jets" 
to distinguish them from the above) are allowed to have an additional relativistic 
velocity (F c 3> 1) relative to the mean jet flow (characterized by Tj). The latter 
situation may be the outcome of relativistic (Petschek-type) magnetic reconnection 
events^ occurring in a highly magnetically-dominated (initial a ^> 1, i.e., magnetic 
energy 3> rest mass energy of the matter) electron-proton jetl 83 ! 84 ! If the fields can 
be efficiently dissipated (downstream er e < 1), with strong mini-jets developing at 
distances of ~ 100r s and reaching r c ~ ^fcr = 10, up-scattering of external nuclear 
photons (e.g., from an ADAF, see Ref . [55) could account for the observed TeV emis- 
sion. Such a model is interesting in that it allows for strong Doppler beaming effects 
and seems capable of accommodating even shorter VHE variability timescaleP^ On 
the other hand, it is not clear whether the assumed conditions are likely to pre- 
vail: Characteristic magnetization parameters for electron-proton (disk-driven) jets, 
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for example, are usually believed to be much smaller (er <C lOOpSl whether/how 
efficient (non-thermal) power-law electron acceleration much beyond the thermal 
Lorentz factor y/am p /2m e ~ 10 4 can be achieved, still remains to be understood. 
Moreover, for the assumed conical jet geometry (opening angle ~ but see 

2.(3) above), the dynamical role of the anticipated guide-field in the jet frame at 
r ~ 100r s is likely to be non- negligible (i.e., the reconnecting fields are not strictly 
anti-parallel), in which case the magnetization of the outflow (<r e > 1) could be still 
high, thereby leading to only weak and slow mini-jets. 

4.3. Hadronic jet models 

Hadronic interactions could in principle also contribute to or even dominate the 
observed VHE emission in M87. 

(1) A hadronic inner jet model including synchrotron emission of energetic primary 
protons as well as the contribution of secondary charged pions (from photo-meson 
/^-production) and muons (from pion decay) has been put forward a while ago to 
describe its time-averaged spectral energy distribution.^ If one assumes proton syn- 
chrotron emission to dominate, maximum (co-moving) proton energies (assuming 
shock- type acceleration with t' acc — T]r' g /c and efficiency rj > 1) would be lim- 
ited to 7^ ~ 4 x 10 10 ?/ _1 / 2 (30 G/B') 1 / 2 . The resultant observed synchrotron emis- 
sion would then essentially be constrained to energies hv p ~ 0. 315/7/ TeV (above 
which the spectrum decays exponentially). Hence, if M87 would be of the mis- 
aligned blazar-type (with small Doppler factor D), one would not expect its hard 
TeV spectrum to be accounted for by proton synchrotron emission. The mean pion 
energy from photo-meson production, on the other hand, is (in the delta func- 
tional formalism) 7^ ~ 0.27^771^/771^, so that hv^ ~ 28hv p ~ 8D/77 TeV. This 
process could thus produce VHE emission up to ~ 10 TeV if protons are able 
to reach maximum Lorentz factors j' p , i.e., provided extreme acceleration 77 ~ 1 
would indeed occur. As r\ oc (c/u s ) 2 , this would require relativistic shock speeds be- 
yond what has been seen so far. It seems not clear yet, whether such a model 
could self-consistently produce the hard TeV spectrum. Photon-meson produc- 
tion in M87 is expected to take place on characteristic (co-moving) time scales 
t' c £ 1.7d {L' 7 /5 x 10 41 erg/s)(5 x 10 15 cm/R') 2 (e 7 /0.01eV), so that the observed 
variability could perhaps be compatible with current constraints. 

(2) Alternatively, VHE 7-rays may originate in inelastic proton-proton collisions 
(with characteristic cooling time scale t' pp ~ /n' p sec). The required high target 
proton density (n' p ) needed to account for day-scale VHE activity could possibly 
be introduced by a star or a massive, dense gas cloud (of typical initial radius 
r c0 ~ 10 13 cm and mass M c ~ 10 29 g) penetrating into the jet on scales of some 
tens of Schwarzschild radii. 88 89 90 If the jet is powerful enough and if a sufficiently 
large fraction (77 > 0.1) of the power tapped by a cloud can be channeled into effi- 
cient acceleration of protons, such a scenario could well explain hard and variable 
VHE emission from non-blazar AGN. For an almost collimated jet (opening angle 
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9^2°) detailed spectral modeling of the M87 VHE flare in April 2010 implies a 
jet power at the upper end of the estimates presented above (cf. 2. (5)) Assuming 
that to be the case, the model can successfully explain the observed VHE (flare) 
light curve and spectrum. To accommodate the observed, large jet opening angle 
(and thus the jet radius Tj, cf. 2.(3)), this approach requires a situation where most 
of the jet energy flux is still concentrated into a narrow core (spine) as otherwise 
the characteristic expansion time scale of the cloud (determining the variability) 
might become too large and the jet power tapped by the cloud (oc (fc/fj)) too 
small to account for the observed VHE characteristics. It seems conceivable that 
such a situation could arise in a hybrid black hole/disk-driven jet scenario. 
(3) A combined lepto-hadronic jet model (including proton synchrotron and proton- 
proton interactions as well as energetic particle transport by convection) has been 
analyzed recently,^ showing that the time-averaged nuclear broadband spectrum 
of M87 might be reasonably accounted for (with its high energy SED part domi- 
nated by pp-processes), if the kinetic jet power would indeed be able to approach 
(perhaps occasionally) values as high as Lj ~ 10 46 erg/s. This seems certainly a 
strong condition when compared with current time-averaged estimates for the jet 
power (cf. 2.(5)). 

A unique signature of a hadronic jet model would be the detection of high energy 
neutrinos. However, typical model predictions lie a few orders of magnitude below 
current instrumental sensitivity. 

4.4. Magneto spheric VHE emission 

The observation that the VHE emission varies on characteristic timescales compa- 
rable to the light travel time across (at most) a few Schwarzschild radii, and the 
finding that the VHE flaring activity in 2008 was accompanied by an increase in 
radio flux close to core, suggest that the relevant non-thermal particle acceleration 
and emission processes could in fact take place in the very vicinity of the central 
black hole itself.^ Magnetospheric models, where these processes are considered 
to occur at the base of a rotating black-hole-jet magnet osphere J3 | 6[l | 1 ° 2 | 104 | 1 03 | ma y 
then become particularly interesting. In principle, at least two main variants can 
be distinguished (see Ref . [3T| for an extended review) : 

(1) Particle acceleration in magnetospheric vacuum gaps (direct electric field ac- 
celeration), with primary electrons in M87 injected by, e.g., annihilation of ADAF 
disk (MeV) photons and quickly reaching high Lorentz factors (j e ~ 10 10 )P^ Cur- 
vature emission and inverse Compton (IC) up-scattering of ADAF soft photons can 
then produce VHE photons with a spectrum extending up to 10 4 TeV. VHE pho- 
tons with energies below some TeV might be able to escape absorption (see below) , 
while IC photons with energies well above 10 TeV are expected to interact with the 
ambient radiation field and initiate pair cascades (e.g., Refs. 11041 11050 just above 
the gap, leading to a large multiplicity.^ The results suggest that this could ensure 
a force-free outflow above the gap and account for the appearance of the VLBA 
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jet. Intermittencies of the cascade process, induced by, e.g., modest changes in ac- 
cretion would give rise to the variability of the TeV emission and to fluctuations of 
the resultant force-free outflow. The HE and VHE parts of the spectrum are then 
likely to be shaped by different emission processes. Models of this type depend on 
a sufficiently spinning black hole to power the VHE emission, require the presence 
of charge-starved regions in the magnetosphere and rely on a radiatively inefficient 
disk configuration to ensure escape of TeV photons. Significant 77-absorption fea- 
tures should become apparent in the VHE spectrum above ~ 10 TeV. 
(2) If the plasma density is such that the parallel electric field becomes effectively 
screened, inertial (centrifugal) effects could probably still allow efficient acceleration 
of electrons to energies of several TeV. The general picture would then be one where 
energetic particles experience the centrifugal force while moving along rotating mag- 
netic field lines anchored in the innermost part of the accretion disk, with electrons 
reaching suitable high Lorentz factors for IC up-scattering of ambient ADAF soft 
photons to the VHE regime close to the light cylinder 77, Models of this type 
still need a sufficiently spinning black hole (Kerr parameter ~ 0.8) for 77 to be 
small enough to account for day-scale VHE activity. In addition, any encountered 
random magnetic fields should be small enough for electron synchrotron emission 
not to exceed current HE constraints! 60 ! 104 ! 

4.5. Escape of VHE j-rays? 

An important issue for inner jet scenarios, and in particular magnetospheric VHE 
models, concerns the question whether TeV photons are able to escape photon- 
photon interaction (77-absorption) in the source's vicinity. In principle, TeV pho- 
tons of energy e interact (producing pairs) most efficiently with target photons 
in the infrared regime ex — (ITeV/e) eV. Escape would thus only be possible if 
the infrared (IR) photon field is weak or sufficiently diluted. This could happen if 
the observed IR field is dominated by e.g. (anisotropic) synchrotron jet emission 
and/or a dusty torus on larger scales. For M87, ground-based (Keck LWS, 300 mas 
resolution) mid- infrared (MIR) observations at 0.1 eV suggest (under the assump- 
tion of isotropy!) a nuclear luminosity Lmir — 10 41 erg/s within 20 pc from the 
black hole (cf. also Gemini/460 mas resolution, and VLT-VISIR/350 mas resolu- 
tion, MIR observations for similar result j£Z| " | 10 o j) IH1 The MIR findings are broadly 
consistent with a non-thermal synchrotron jet origin of the nuclear IR emission.^ 
Along with other, more recent IR observations (Spitzer MIR 1.5 — 6 x 10 13 Hz, 
~ 3" resolution; Herschel FIR 0.6 — 3 x 10 12 Hz, ~ 8" resolution), they seem to 
disfavor a significant /dominant (quasi- isotropic) IR contribution by thermal dust 
(torus) emission in M87 (but see also Ref.[SE]fb r a possible dust contribution towards 
longer wavelengths) Experimentally, this would imply an optical depth for 

TeV photons r 77 (e) ~ a^n s {e T )R d = 0.0577( J L T /10 41 erg/s)(0.1pc/R d )(e/lTeV), 
considering the IR contribution to be dominated by jet emission (with associated 
anisotropy rj <C 1) on scales Rd ~ 0.1 pc. Under such conditions TeV photons may 
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well escape absorption. Alternatively, one may calculate the expected IR photon 
field in the vicinity of the black hole by specifying some inner disk configuration 
and accretion rate, and then use this to derive optical depths. If the inner disk would 
be of the standard (i.e., cooling-dominated, radiatively efficient) type, the emergent 
disk photon field would peak in the IR and be dominated by emission from regions 
close to the marginal stable orbit, thus making an escape of TeV photons from 
regions close to the black hole rather unlikely! 31 ! 101 ! Most probably, however, the 
(inner) accretion flow in M87 is radiatively inefficient and dominated by advection 
(ADAF, see above 2.(4)). In such a case, escape of TeV 7-rays would again be 
feasibleP^ Even if, for example, the whole nuclear radio to X-ray flux observed in 
M87 would solely originate in an ADAF (which is certainly over-restrictive given 
that most of its nuclear radio-X-ray emission is usually believed to be jet-dominated, 
e.g. Ref. l97|) . then in the presence of a rapidly rotating back hole ~ 10 TeV pho- 
tons would be able to escape absorption once produced on scales ~ 5r s P^l Detailed 
calculations show that the anticipated optical depth for TeV photons is sensitively 
dependent on accretion rate) 30 l 93 J so that magnetospheric VHE emission (provided 
it indeed occurs) may in fact be modulated by small changes in accretion conditions. 
Given current evidence there are thus good reasons to believe that, at least for M87, 
TeV 7-rays may be able to escape from the vicinity of its supermassive black hole. 
This in turn suggests that future high-sensitivity instruments such as the European 
CTA project could play a particularly important role in testing magnetospheric 
VHE 7-ray emission scenarios in M87. 

4.6. On dark matter, ultra- high- energy cosmic rays and EBL 
prospects 

The observed VHE variability excludes a dominant gamma-ray contribution due 
to cosmic-ray interactions with the thermal gas in the interstellar mediurrP^fil or 
neutralino annihilation inside the dark matter halo of M87, 107 although it remains 
possible that some of these processes could contribute to its low state (long-term 
average) gamma-ray emissionP^I 

The relevance of M87 as candidate source for the origin of ultra-high-energy cos- 
mic rays (UHECR, > 4 x 10 19 eV) has often been discussed in the literature, see 
e.g. Ref. I109[ IUCH 1111! 11121 Efficient proton acceleration to ultra-high energies in 
magnetospheric vacuum gaps, however, appears rather unlikely (due to curvature 
losses and partial field screening) | 30 l 31 l On the other hand, diffusive shock accel- 
eration of UHE protons in the jet of M87 (cf. also the proton synchrotron model 
described above) seems marginally possible provided relativistic shock speeds would 
indeed occur. Given its proximity and anticipated relatively small deflection angles 
(for protons), M87 is one of a few potential extragalactic sources that could possi- 
bly allow to test UHE cosmic-ray acceleration by directional ~ 10 19 eV cosmic-ray 
astronomyl 114 ! 115 ! 

The H.E.S.S. observations in 2005 indicate that the VHE spectrum of M87 could 
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well reach up to a few tens of TeV. Due to its proximity, the VHE 7-ray flux of M87 
is likely not to become significantly affected by interactions (pair production) of TeV 
photons with infrared photons of the diffuse extragalactic background light (EBL), 
at least up to some tens of TeV. An extension of the hard power-law 7-ray spectrum 
of M87 beyond 30 TeV, for example, would thus allow a very unique probe of the 
EBL in the transition region from mid- to far infrared (A ~ 30 [e 7 /30 TeV]/zm)P"^ 
The measured TeV fluxes and anticipated sensitivities of the future Cherenkov array 
CTA could make this a prime physics case. 

5. Conclusions 

The radio galaxy M87 is a unique source on the extragalactic sky. Its proximity 
and huge central black hole mass allow for an analysis of fundamental astrophysi- 
cal processes in unprecedented detail. M87 is well-known for a one-sided, kpc-size 
relativistic jet visible at radio, optical and X-ray wavelengths. However, compared 
to typical blazar sources, this jet is most likely substantially inclined ~ (15 — 25)° 
with respect to the observer, so that the broadband emission observed from M87 
is not necessarily (everywhere) swamped by relativistically-beamed jet radiation. 
M87 thus appears to be a special active galaxy where we may be able to discover 
central features related to other high-energy processes. The rapid (days-scale) VHE 
flux variability detected during high source states by current 7-ray instruments is 
the fastest variability seen so far at any wavelength, and suggests that the variable, 
hard-spectrum 7-ray emission may originate close to its supermassive black hole. A 
variety of physical scenarios have been proposed in the literature to accommodate 
such findings. Some of the pros and cons/limitations of these candidate models are 
mentioned in this review, and the relevance of further VHE observations is stressed. 
VHE flaring characteristics (energy spectra and extension, variability time scales) 
like those reported bear the rare potential to give us important insights into the 
near black hole environment of active galaxies. Future, more sensitive VHE ob- 
servations (CTA) will allow to study spectral variability and further constrain the 
time scales of VHE flux variations - both important inputs for theoretical model- 
ing. Simultaneous VHE/high- resolution radio observations, on the other hand, will 
have the potential to pin down the location of the VHE emission on spatial scales 
~ 10 2 r,,.. M87 is the best-suited AGN for such kind of studies, since it is (beyond 
Sgr A*) the most promising object for which direct imprints of the black hole in 
future high-resolution radio observations can be expected.^ 
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Fig. 1. The large-scale jet of M87 as seen in radio, optical and X-rays. The bright knot A is 
located at about 1 kpc from the nucleus (Credit: X-ray: NASA/CXC/MIT/H. Marshall et al., 
Radio: F.Zhou, F.Owen (NRAO), J.Biretta (STScI), Optical: NASA/STScI/UMBC/E.Perlman et 
al.) 
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Fig. 2. The small (sub-parsec) scale jet in M87 as observed at 43 GHz with VLBA in 2007. The 
resolution achieved is 0.43 by 0.21 mas. Note that for the distance (16.7 Mpc) of M87, 1 mas = 
0.081 pc, and 100 Schwarzschild radii=0.36 mas (for M B H = 3 X 10 9 M Q ). (Credit: NRAO: R.C. 
Walker et al. , |http:/ /www.aoc.nrao.edu/~cwalker/M87[ l 
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Fig. 3. VHE flare of M87 as observed in April 2010 by different instruments. Day-scale activity 
is evident. The curve shows a fit of an exponential function to the data. From Ref. ^) 
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Fig. 4. Characteristic (non-simultaneous) nuclear SED of M87. A simple HE (Fermi) power-law 
extrapolation could fit the VHE (H.E.S.S.) low state data, but would under-predict the VHE high 
state data. Historical data are, in increasing order: Radio and sub-mm data [green, filled circles]: 

5 GHz VLBI (0.7 mas resolution, Pauliny-Toth et al. 1981), VLBA (sub-mas resolution, Zavala 

6 Taylor 2004), VLBI 22GHz (0.2 mas resolution, Spencer & Junor 1986), IRAM/89 GHz (2 
arcsec resolution only, upper limit for nuclear flux, Despringre et al. 1996), VLBI/100 GHz (0.1 
mas resolution, Baath et al. 1992), SMA 230 GHz (300 mas resolution, Tan et al. 2008). Spitzer 
FIR data [green, open circles] (several arcsec resolution, Shi et al. 2007). MIR data (cyan, star) 
from Keck (~ 300 mas resolution, Whysong & Antonucci 2004) and Gemini (460 mas resolution, 
Perlman et al. 2001). Optical HST data (pink, squares) are from Sparks et al. 1996 (~ 200 mas 
resolution). Nuclear X-ray data (black, diamonds) are from Chandra observations (~ 500 mas 
resolution, Marshall et al. 2002) in July 2000. Note, however, that the nuclear X-ray flux may 
increase by a factor ~ (5 — 10) during strong flaring states (as e.g. in February 2008, Harris et 
al. 2009). Swift/BAT hard X-ray longterm (03/2005-01/2009) upper limits [grey, open triangles] 
(Ajello et al. 2009; Abdo et al. 2009). Fermi-LAT longterm HE data (08/2008-05/2009) [orange, 
filled squares] (Abdo et al. 2009). VHE data from H.E.S.S. for the 2005 flaring (blue, triangles) 
and 2004 low state (red, filled circles) (Aharonian et al. 2006).) 



